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I. INTRODUCTION 

Crucial features of QCD are believed to be encoded in 
the infrared behaviour of its Greens functions. Indeed, 
the greatest unsolved problems of the theory - confine- 
ment, dynamical chiral symmetry breaking and the emer- 
gence of a mass gap - presumably have their origin in 
the infrared, where interactions are strong, the coupling 
is large and perturbation theory breaks down. 

Thus it is extremely valuable to have methods at hand 
which allow us to study this sector. As argued in [1, 2], in 
the deep infrared, far away from all other scales (Aqcd 
or quark masses), the theory should be conformal and 
propagators D should thus exhibit an asymptotic power- 
law behaviour, 



2\<S t 



Di(p) ~ (p 2 ) 



(1) 



with infrared critical exponents <5; (and coefficients Cj 
which contain some power of a renormalization scale 
//.) These exponents can be extracted from the Dyson- 
Schwinger equations (DSEs), the equations of motion for 
a quantum field theory. 

This endeavor has been pursued with great success in 
the Landau gauge [3], where, in accordance with Gri- 
bov's confinement scenario [4-6], the infrared suppres- 
sion of the gluon propagator could be traced back to the 
enhanced divergence of the ghost propagator. 1 

Unfortunately, other gauges seem to be more difficult 
to access. In the Coulomb gauge, which is of particular 
interest for the present authors, the SO(3,l) symmetry of 
spacetime is reduced to SO(3), thus a propagator D(p) 



generically depends on two physical momentum scales p\ 
and p 2 . 

Accordingly various different infrared limits have to be 
distinguished: Even if one has p 2 — > and p% — > 0, one 
could have completely different behaviour depending on 

the ratio py. 2 Even in a quasi-instantaneous approxi- 
mation (where all diagrams without at least one instan- 
taneous propagator have been neglected) infrared criti- 
cal exponents have turned out to be rather elusive ob- 
jects [10]. 

Thus in this article, instead of working with the 
ground-state theory, we study QCD (or rather SU(AT) 
gauge theory) in the limit of very large temperature, 
where certain simplifications occur. At first glance one 
may wonder how useful this could be, taking into account 
that according to common lore, at large temperatures 
QCD becomes "deconfined" , chiral symmetry is restored 
and the mass gap is gone. 

However, as discussed in section 1 of [11] (where addi- 
tional references are given), there is increasing evidence 
that this picture is likely to be incomplete or even wrong. 
Perturbative [12, 13], lattice [14], functional [15, 16] and 
experimental [17] results suggest that the infrared sector 
and bound states play an essential role also at very high 
temperature, so that there is no "deconfined phase" in 
the strict sense. 

Thus one can hope to obtain useful information about 
the infrared sector even in the case of very large temper- 
atures, which we will study in the following. 



Organization of this Article 



1 Note that some recent lattice studies question these results. 
While those concerns are certainly to be taken seriously, it would 
seem overly hasty to dismiss the results obtained by functional 
methods. In particular, since the finite size of a lattice corre- 
sponds to an infrared cutoff in momentum space, huge lattices 
are required to perform a reliable extrapolation to infinite vol- 
ume, and there is still an ongoing debate about systematic er- 
rors [7, 8]. In general, there is a discussion about the role of the 
"massive" as opposed to the "scaling" solution [9] . 



The article is organized as follows: In sec. II we give 
the action of Coulomb-gauge Yang-Mills theory for fi- 



2 In general one encounters all problems present for func- 
tions of two real variables, where not only in general 
lim^o lim H ^o f(x, y) ^ lim^-,0 lima;->o f(x, y), but the limit 



Urn 



f(x n , y n ) may even depend on the precise path 



on which the origin is approached. 
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nite (and in particular extremely high) temperature; in 
sec. Ill we define propagators and proper two-point func- 
tions and discuss the relation between them. 

In sec. IV we state the Dyson-Schwinger equations 
of the theory, which are examined more closely in a 
bare-vertex truncation. Since these equations contain 
infrared-divergent integrals, we discuss the topic of these 
singularities in sec. V. 

We define the infrared critical exponents in sec. VI and 
employ power-law ansatze in sec. VII in order to obtain 
the infrared asymptotic Dyson-Schwinger equations. 

The exponents are constrained analytically and deter- 
mined numerically in sec. VIII. (The evaluation of power- 
law integrals is discussed in appendix A.) In addition we 
can also determine relations between the propagator co- 
efficients, which are derived in sec. IX. 

In sec. X we try to estimate the range of validity of the 
approach followed in this article, and in sec. XI we discuss 
the results obtained for the infrared critical exponents 
and give a brief summary. 



II. LOCAL ACTION 



We rewrite the theory in the first-order formalism, 
by introducing a new field 7r" by a Gaussian iden- 
tity, so the action reads 



1 

T 



d s x 



mi(-DiAo) + -7T? + -F§ + dicD lC 



(5) 



The new field can be interpreted as the momentum 
conjugate to Af and thus plays the role of a color- 
electric field. It can be decomposed into transverse 
and longitudinal parts 



7T- - di<f, 



(6) 



where 7r£ is transverse, = 0, which gives the 
action that will be used to derive the DSEs, 



b ~ T 



d s x 



ifa-M(-A4>) + 5fa) a 



(7) 



Our starting point is the Yang-Mills action in d = s + 1 
dimensions, gauge-fixed to the Coulomb gauge, which is 
particularly well-suited for studies of finite-temperature 
field theory, 



S 



YM.Coul 



= d S " 



(2) 

where F M „ = d^A v - duA^+gA^ x A v , and (A^ x A v ) a = 
f abc A^AI; the gauge-covariant derivative is given by 
Die = Di[A]c = die + gAi x c. We will now modify 
the action (2) in four ways: 

• We apply the on-shell formalism, so the Nakanishi- 
Lautrup field b is integrated out in order to directly 
impose the transversility condition 



diAi = . 



(3) 



We turn to finite temperature, so the temporal in- 
tegral has the limits J Q dx , where /3 = ^ and 
T is the temperature. Integrals over fco will be 
replaced by a sum over Matsubara frequencies, 
/d**>-T£ n . 

We neglect all but the 0th Matsubara frequency, 3 
which is the same as dropping all time derivatives, 

8q — > 0, in the action, and replacing f dxo — > 
so the action simplifies to 



Sl = ? / dS ' T (l {DlAo)2 + 1*% + 9fcD ^ 



(4) 



3 Note that the Linde problem [12] has its origin in the zeroth 
Matsubara frequency as well. 



III. DEFINITION OF PROPAGATORS AND 
PROPER 2-POINT FUNCTIONS 

If we confine ourselves to the zero Matsubara fre- 
quency, propagators only depend on the spatial momen- 
tum. The propagators of the transverse fields are defined 

as 



(A«{x)A)(y)) 
«(zK-0/)> 



d s k 
d s k 



e ik < x -y^6 ab pT(k)D AA (k), 
jH*-v) 6 <*pT(k)D wn (k), 



where P^(k) is the transverse projector, 



(8) 



(9) 



The propagator (w' i (x)Aj (y)) is proportional to fco both 
at tree-level and for the power-law ansatze employed in 
sec. VI. Thus it has vanishing zero-Matsubara compo- 
nent in this context, and will be neglected in the asymp- 
totic infrared limit, 



(n' i (x)A j (y)) = 0. 



(10) 



This removes the mixing of the transverse fields, so the 
proper functions are given as the one-dimensional inverse 
of the propagators, 



r AA (fc) 



i 



£>aaAt'*' - D 2 A ^ D AA {k) 
Daa 1 



DaaD^-D a ^ IW(fc) 



(11) 
(12) 
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On the other hand the scalar Bosc fields do mix. Their 
propagators are defined by 



(A a (x)A b (y)) 



(A a (x)<p b (y)) = 



d s k 

d s k 

d s k 



jH*-v) S <*D AoAo (k), 
e ik < x -^S ab D Aoip (k), 
e ik -( x -y)S ab D V(p (k) , (13) 



the Faddeev-Popov ghost propagator is defined by 

d s k 



(c a (x)c b (y)) 



(27T) 



e M^-v)s ab D C c{k). (14) 



While the inversion of the ghost propagator (in order to 
obtain the proper 2-point function) is simple, 



r Ec (*0 



i 



for the other scalar fields the proper 2-point functions are 
two-dimensional matrix inverses of the propagators, 



r A a A (k) 



D VV ( k ) r _ D A A (k) 



^A ov >( k ) = - 



r w (fc) — 

D<pA (k) 
A(k) : 



where 



A = D vv Da a - D 



A ip- 



(16) 



(17) 



IV. 



TRUNCATED DYSON-SCHWINGER 
EQUATIONS 



The derivation of the Dyson-Schwinger equations 
(DSEs) for the theory described by (7) is straightforward, 
yet tedious. (Details of the derivation will be discussed 
in [18].) To simplify this endeavour we neglect the cubic 
and quartic pieces of Fij in (7), since the scalar fields Aq, 
(p, c and c are expected to be dominant in the infrared; 
loops containing a three- or four-gluon vertex (and the 
the corresponding amount of transverse propagators) are 
supposed to be subdominant.. 

To further simplify the equations (and since little is 
known about the dressed vertices of this theory anyway), 
we employ a truncation in which dressed vertices are re- 
placed by bare ones. (Note that in [10] a more general 
ansatz for the vertices did not change the general pic- 
ture.) A graphical representation of the resulting DSEs 
is given in Figure 1. 

With these truncations, the equation for Taa reads 

P?(k)T AA (k) = P?(k) k 2 +g 2 TNPl l {k) I mn P?(k). 

(18) 



where the temperature T stems from the Matsubara sum 
(of which we keep only the zeroth term). The loop inte- 
gral I m n(k) is sandwiched between transverse projectors 
P T {k) 1 and is given by 

d s p 



(2tt)' 



Ar m 7r„ (P) Da A (j> + k) + 

PmPn 



x ^D ipAo {p)D AoV {P + k) + D c - C {p) D c£ (p + k) 

(19) 

where the propagator of the color-electric field has the 
decomposition 

Ar m 7r„(p) = Pl n {p)D^«,{p)+p m p n D w {p). (20) 

The DS equation for reads 

pT(k)T v , v ,{k) =pT(k)+g 2 TNP? m (k)J mn pT(k) (21) 

where the loop integral is given by 



d s p 
(271-)* 



P? nn (p)D AA (p) D AoAo (P + k). (22) 



We further obtain the equation for Ta a 
T Ao A {k)=g 2 TN 



PT(p)D AA (p)D WjVi (p + k), 
(23) 

where D 7T 7ri is given in (20). The DS equation for T w 
reads 

d s p 



r w (fc) = k + g TN 



(27T) 

x D Ao a (p + k) 
the DS equation for T v a is given by 



kikjPlMD AA {p) 



(24) 



iWfc) =ik 2 +g 2 TN 



PT(p)D AA (p) 



(2tt) s ^ y 
xkik J D Aol p{p + k). (25) 
Finally the equation for r gc reads 

d s p 



T 5c (k) 



g TN 



(2n) s 



PMp)D AA (p) 



xhkjD^ip + k). (26) 

The tree-level terms in equations (18) to (26) do not 
directly affect the infrared asymptotic behaviour for es- 
entially two reasons: 

• The scalar fields are expected to be infrared- 
enhanced. Accordingly we impose the horizon con- 
dition [5, 6] on the Faddeev-Popov ghosts and 
(since we expect at least qualitatively analogous 
behaviour) also on the bosonic fields. As a conse- 
quence, the tree-level part in (23) to (26) is can- 
celled by quantum fluctuations. 

• The DS equations for the transverse fields contain 
at least one (uncancelled) loop with at least one 
scalar propagator, which will - due to infrared en- 
hancement - dominate the tree- level part. 




FIG. 1: The system of bare vertex truncated Dyson-Schwinger equations in the first-order formalism, where the tree-level terms 
have (for the transverse fields) been neglected or (for the scalar fields) removed by imposing the horizon condition. Prefactors 
and signs have been absorbed in the graphs. Diagrams and equations drawn in gray drop out in the approximation (10). 
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V. INFRARED DIVERGENCE IN DSE 

In Coulomb gauge, the time-time component of the 
gluon propagator has an instantaneous part, 

-g 2 D AoAo (x, t) = V cou i(r)6(t) + V non 

—inst. t), (27) 

known as the color-Coulomb potential V cou \(r), where 
r = \x\. When the gauge-invariant potential introduced 
by Wilson is confining, linv^oo Vw(r) = oo, it provides 
a lower bound on the color-Coulomb potential asymp- 
totically at large r, 14 ui(f) > Vw(r), summarized by 
"no confinement without Coulomb confinement" [19]. 
In the confining phase we expect the Wilson potential 
to rise linearly, Vw(r) ~ or, where a is the physical 
string tension, and in this case the color-Coulomb po- 
tential rises (at least) linearly, V QQU \(r) ~ ir cou \r, where 
Ccoul is the Coulomb string tension and er cou i > a. 
Moreover in Coulomb gauge gAo is a rcnormalization- 
group invariant [19], which implies that V cou i(r) is also a 
renormalization-group invariant. Thus, as long as cr cou i is 
finite, it has a well-defined physical value in MeV which 
has been calculated in lattice gauge theory [20-24] . 

The linearly rising potential r corresponds in momen- 
tum space to 1/fc 4 , so one encounters infrared divergences 
of the type J d 3 k/k 4 in the DS equations. To address this 
problem, wc observe that the DS equations arc originally 
derived in position space [from the functional identity for 
ST/SA(x)], and moreover these equations remain free of 
infrared divergences when written in position space, even 
in the presence of long-range potentials. Indeed a loop 
integral such as 



L(p) = 



d s k 
(2tt) s 



D 1 {p-k)D 2 (k), 



(28) 



which is a convolution in momentum space, corresponds 
in position space to the ordinary product 



L(r) = £>i(r)£» 2 (r). 



(29) 



This product is well-defined for long-range potentials 
such as Di(r) ~ r. This is the basic observation which 
gives a well-defined meaning to the DS equations for long- 
range forces. 

One has the option of working entirely in position 
space. However, as a matter of convenience, we may 
also work in momentum space, as usual, once we provide 
a well-defined two-way translation between position and 
momentum space. 

The Fourier transform of a long-range correlator such 
as Di(r) = cr cou \r is well-defined by providing a conver- 
gence factor, for example exp(— er 2 ). Indeed the integral 



D(k) 



d 3 x <j cou \r exp(— ik ■ x) exp(— er 2 ), (30) 



has a finite limit for e — > 0, 

D(k) = -87r<7 coul /fc 4 . 



(31) 



The minus sign here is as it should be because cr C oui^ is 
an attractive potential. [The minus sign in (27) was in 
fact introduced to make V CO ui(r) a positive quantity.] We 
conclude that there is no difficulty in taking the Fourier 
transform of long-range potentials. 

However the inverse Fourier transform poses an appar- 
ent difficulty because it has an infrared divergence of the 
form J d 3 k/k 4 . Moreover the DS equations which con- 
tain a loop integral (28) have infrared divergences of this 
type. 

To address these problems, consider the inverse Fourier 
transform in s spatial dimensions 



d s k exp(ik-x) T(s/2-a) 



„2q- 



{2tt) s {k 2 Y 



2 2a W 2 T(a) 



(32) 



It is free of infrared divergence provided that the param- 
eter a satisfies the bound a < s/2. The integral has been 
evaluated by Gaussian integration after insertion of the 
identity, 



1 



(k 2 



T(a) 



dp /T" 1 exp(-/3£; 2 ). 



(33) 



[By this method of integration the convergence factor 
cxp(— ek 2 ) is not needed explicitly] The last integral is 
convergent only for a > 0. However the the left hand side 
of (32) is well defined for all real a satisfying the bound 
a < s/2 so, by analytic continuation, eq. (32) holds for 
all a that satisfies this bound. Because \/ r s ~ 2a appears 
on the right hand side of (32), the restriction a < s/2 
corresponds to a negative power of r. 

To extend the inverse Fourier transform to longer range 
potentials such as r itself, which corresponds to a more 
infrared-singular power of k, we observe that such a more 
infrared-singular power may be written as a derivative, 



1 



-1 



1 



(k 2 ) 



2U+Q 



4ot(*/2 - a - 1) dkf (k 2 ) 



2\a 



(34) 



and we take the inverse Fourier transform of the right 
hand side in the distribution sense. This gives 



d s k 



exp(ik ■ x) 



1 



-1 



d s k 



4a(s/2 - a - 1) J (2tt) s 
i- d s k 



d 2 

cxp(ifc • x) -j^- 2 



1 



4a(s/2 - a - 1) J (2n) 



exp(ifc • x) 



dk'f {k 2 ) a 
1 



2\a 



r(s/2-a-l) 2(a + l)-.s 



2 2(«+l) n s/2 p( Q + X ) 



(35) 



where we have used (32) for a < s/2, and r = (a; 2 ) 1 / 2 . 
We now observe that (35) has the same form as (32), with 
the substitution a — > a + 1. Thus the result of giving 
l/(fc 2 )" +1 a meaning as a distribution is quite simple: 
formula (32) is continued analytically from a to a + 1, so 
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it is valid under the weaker restriction a < s/2 + 1. By 
induction, formula (32) may be continued to a < s/2 + n, 
where n is any integer. 

Finally we note that the loop integral (28) that ap- 
pears in the DS equation becomes well-defined by the 
same method. Indeed, suppose that Di(fc) = l/(k 2 ) a+1 
is written as the derivative 



= E x (k) 



(36) 



where E\(k) is written above and is less singular than 
Di{k). Then the ill-defined loop integral (28) may be 
replaced by the well-defined expression 



m = 



d 2 



d s k 



m J (2tt)< 



E 1 (p-k)D 2 (k). 



(37) 



However there is no need to do this explicitly because the 
result of this substitution may be obtained, as before, by 
analytic continuation in a of the original loop integral 
J^Jj-D 1 (p-k)D a (k). 

We conclude that instead of working in position space, 
where there are no infrared divergences, we may, as a 
matter of convenience, work directly in momentum space 
according to the following prescription: The standard 
loop integral (28) is performed for values of the critical 
exponents a for which it is well defined. The result is 
then analytically continued in a to the values of interest. 



VI. DEFINITION OF INFRARED CRITICAL 
EXPONENTS 

The only dimensionful parameter in the DSEs is g 2 T, 
which in spatial dimension s provides a mass scale m 
defined by 

4— s 2rr 

m = g T. 



(38) 

As an ansatz we look for a solution to the DSE for 
the one-Matsubara frequency propagator that is a simple 
power law in the spatial momentum, 



b*m 2aA 



2U + QA 



ZW(fc) 



{k 2 ) a «> 



b gll m 2a ^ b Q m 2a ° 
D cc(k) ~ „„ 2 m+„„„ ; D AoAo (k) ~ „^,+ nn ; (39) 



D <fiA (k) 



(fc2)l+a gh 

-ib m m 2a ™ 
(k 2 ) 



2\l+a„ 



(k 2 y+ a °' 

h in ^ a ^ 



(k 



In the following it will be sometimes convenient to write 
av = I + Qtt' ■ The mass m cancels out of all equations 
because of engineering dimensions. 



to the other scales in the theory, and we take the in- 
frared asymptotic form of the propagators. This will 
yield a finite system of equations which will provide a 
self-consistent infrared limit of the propagators. 

The infrared asymptotic equations read (with loop in- 
tegrals Is(v,s)-, Is(v,v) etc- to be defined in subsec- 
tion VII B) for T^^ 4 

u 

b A b<pI S (v,S)(otA, a v ) 



bob v + b 2 m 



for r 



<pA ' 



bob v + b 2 



+ &a67t'/s(V,V)( q A, av)i (40) 



= —bAb m Is(y,S)(<XA,a m y, (41) 



for r 



tptp, 



, 7 b °, 7 2 = b A b I s{v ^ s) (a A ,a ); (42) 

for T A a, 

1 = b A bobn>Iv(v,s)(&iv',ao) + bAbob ip Iv(s,s)(ao-,a<p) 
-b A b m I V (s,S){oi m ,a, m ) + b A b sh I v ^s,s){otsti, a g h)- 

(43) 

for IV ^ , 

1 = b A bob v 'I V (y^s){a A , a ); (44) 

for r & , 

1 = -b A bl h I S (y,s)(oiA, a g h) ■ (45) 

A. Symmetry of the infrared asymptotic equations 

There is a two-parameter continuous symmetry trans- 
formation that these equations inherit from the cubic in- 
teraction terms irigAi x Aq and dicgAi x c, namely 

A l — > exp{i(3)Ai; A — > cxp(i7)A ; c — > cxp(i7)c 
u l -> exp[-i(/3 + 7)]7r;; c -> cxp[-i(/3 + j)]c. (46) 

As a consequence of this symmetry, the infrared asymp- 
totic DS equations are invariant under the transforma- 
tions of the asymptotic propagators 

b A -> exp(2i/3)6 A ; K> -> exp[-2i(/3 + 7)] V 

6 g h — > cxp(-i/3)&g h ; feo — > exp(2i7) 6 ; (47) 



exp(-i/3)6 m ; 



exp[-2i(/3 + 7 )] 6, 



VII. INFRARED ASYMPTOTIC DS 
EQUATIONS 

In the DS equations, we take the external momen- 
tum k and the loop momentum p to be small compared 



4 We have checked explicitly that in the infrared, given the DS 
equations of Figure 1, both the assumption D vv Da a q •C ^\ ov 
and D vv Da Aq 2> ^\ ov lead to a contradiction for this system 
of equations. Thus D vv Da a an d ^\ ov have the same in- 
frared behavior. A possibility we have not further explored in 
this article is a cancellation bo?V + = 0- 
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Because of this 2-paramcter symmetry the DS equations 
provide only 4 relations among the 6 6-cocfficicnts. 



B. Definition of loop integrals 

We now define the symbols that represent the loop 
integrals, 

I s{v , s) (a A ,a v ) = iV fc - s + 2 «A+2a v+ 2 (48) 

d s p k 2 p 2 — (p ■ k) 2 



where 



Iv(v,S)(a A ,a ) 



(27r) s (p2)2+a A [(fc-pjajl+a* 
J\TU— s+2«A+2ao+2 



S - 1 

d s p (s - 2)k 2 p 2 + (p- k) 2 



(49) 



^s(v,v)(a A ,aV) = Nk 



{2ir) s (p2)2+a A [(fc_p)2]l+ a o 
s+2q a +2q,,+4 / d s p .^p. 



(2tt) s 

( S ~2)(fc-p)V + [(fc-p)'P] 2 



( p 2)2+ QA [(/c_p)2]2+a x , 

_/yi,-s+2Qo+2a ¥ , 
s - 1 

ePp fc 2 p 2 — (p- k) 2 



(51) 



(2tt) s (p2)l+a p_ p )2]l+a v 



One sees by inspection that two of the symbols are related 
by 

I v(s,S){^iOt v ) = (s - l) _1 / s( y iS) (Q;o - l,a v ). (52) 

The symbols have the value (see appendix A for the 
exemplary evaluation of one such integral) 

N(s-l) T(2 + a A + a gh -s/2) 
I S( vM*A,a gh ) = r(2 + aA ) 

r(s/2 - a gh ) T(.s/2 - a A - 1) 



T(l + agh) T(s - a gh - a A - 1) ' 

(53) 

T , , Net T(2 + a A + a -s/2) 
/ ms) (a A) ao) = (54) 

;; r(s/2-a -l) r(s/2-a A -l) 
r(l + ao) r(s — an — a A — 1) 

where 

ci = (s-l)(4a A + 3)-(2a A + 2ao + 3)(2a A + l), (55) 



^s(v,y)(aA,a 7r /) 



Nc 2 r(2 + a A + 6y - a/2) 
4(4tt) s /2 r(2 + a A ) 

T{s/2 - ay - 1) r(s/2 - a A - 1) 
r(2 + av) r(s-a A -ay-2) 

(56) 



c 2 = (s - 1) [s - 1 + (1 + 2a A )(l + 2av)] (57) 



4-(s,s)(ao,av) 



AT r(l + a + a v - s/2) 



2(4tt) s / 2 T(l + a ) 

r(s/2-a y ) r(s/2-a ) 
r(l + a v ) r(s — a^, — ao) 



.(58) 



C. Check of loop integrals 

We obtain a useful check on the evaluation of the loop 
integrals by rewriting 5 the numerator of the integrand of 



I 



s(y,v). 



Af = Ti-[(q 2 -qq)(p 2 -pp)}, 
where q — k — p. We have 

T = q 2 (.s-l)p 2 -q-(p 2 -pp)-q 



(59) 



a 2 

1 1.2 



k 2 



k (s — l)p — k ■ (p — pp) ■ k 



= |j Tr[(fc 2 - kk)(p 2 - pp)] 



k 2 



- 1 ) k ■ (p -pp)-k, 



(60) 



which, by comparison with (48) through (51), leads to 
the identity 

Is(y,V){<XA,6iir') = (s - l)I V ( V)S )(a A , 

+Is(V,S)(<XA,&ir') (61) 

-(s - l)/y ( 5 iS) (l + a A , 1 + av). 

As a precise check, it has been verified that this relation 
between the 4 integrals is satisfied by the 4 values just 
given. 



VIII. DETERMINATION OF INFRARED 
CRITICAL EXPONENTS 

A. 4 power-relations among infrared critical 
exponents 

There are 6 infrared critical exponents and 6 DS equa- 
tions. By equating powers of momentum on both sides 



5 In the following the expression p\p2 denotes the dyadic product 
of the s- vectors pi and p2- In a less compact way, one would 
write (59) as 

M = Tr s [(|q| 2 l s - <jq T ) (|p| 2 l s - pp T )] , 

where p and q are column vectors and T denotes transposition. 
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of the DS equations, one obtains relations between the 
infrared critical exponents. From the equation for r ec , 
one obtains 

2a gh + a A = (s - 4)/2; (62) 

from the equation for T^, 

a + a v + a A = (s - 4)/2; (63) 

from the equation for T v a , 

2a m + a A = {s- 4)/2; (64) 

from the equation for Y^'ir', 

Q'o + QV + a A = (s - 4)/2. (65) 

These 4 power relations come from the 4 DS equations 
that have only one term on the right-hand side. They 
leave undetermined two infrared critical exponents which 
we may choose to be a g h and ao • The remaining infrared 
critical exponents may be expressed in terms of these by 



The equation 



a A = (s - 4)/2 - 2a gh 
o-k' = a v = 2a g h — ao . 

a + a v = 2a m , 



(66) 



(67) 



which follows from the above, relates the critical ex- 
ponents of D AoAa ~ l/k 2+2a ", D w ~ \/k 2+2a * and 
D V a ~ l/k 2+a f +a °. Thus the infrared critical expo- 
nents «o an d ct v characterize the elementary fields A 
and (p. 



This equation allows us to write (40) as 
b v = b A b v b 2 gh I s{ y. s) (a A ,a ip ) + b A b 7Tl b 2 ^I s{y y ) (a A ,a^) 



or, by (45), 

V _ -^s(v,s)(«A,a g h) - Is(v,s){a^a v ) _ 
b v Is{v.v){a A ,a^) 



(70) 
Fi. (71) 



Likewise (69) allows us to write (43) as 
1 = bAbob v >I V (v,s)(&n',aio) + bA.bob S pI V (s,s)(ao,a< p ) 



-b A (bob v - b sh )Iy(s,S)(°tm, Oi m ) 
-&A& g h^y(S,S)(«gh, CUgh)- 



(72) 



With a m = a g h, there is a partial cancellation between 
the last two terms which are the contribution from bose 
and fcrmi ghost loops respectively, and we obtain 

1 = bj^bobn'Iy^y s^an', a ) (73) 
+b A bob v [I v{SiS - ) (a ,a v ) + I V ( S ,S)(<Xgk,<Xsb)]- 

This gives, by (44) 

V _ Iv(s,S)(oto,a<p) +-?y(S,s)(agh,a g h) _ 



\ I V (y S )(a A ,a ) - I v(y ^(a n ',a ) 



F 2 . (74) 



From equations (71) and (74) we obtain the final equation 
that determines the infrared critical exponents. 



F x -F 2 = 0. 



(75) 



Equations (68) and (75) together with the above power 
relations given previously determine the remaining two 
infrared critical exponents, ap and a g h- 



B. Equations for ao and a 



C. General Remarks on the Equations 



When the above 4 power relations on the 6 critical 
exponents (the 6 a's) are satisfied, the power relations 
among the remaining 2 DS equations are satisfied identi- 
cally. The 6 DS equations also provide 6 relations among 
the 6 6-coefficients. However, because of the 2-parameter 
symmetry invariance (47), of these 6 equations, only 4 
are independent conditions on the 6-coefficients. The 
remaining two equations provide consistency conditions 
that determine the two missing relations among the in- 
frared critical exponents, as we now show. Thus all 6 
infrared critical exponents are determined. 

From (41) and (42) we obtain an equation relating crit- 
ical exponents, 



Is(v,s)i a A,at gh ) — —I S ( ViS )(a A ,a ), 



(68) 



where we have used a g h = a m . 

The remaining relation between critical exponents is 
obtained as follows. From (41) and (45) and power rela- 
tion a g h = a m we obtain 



62b 



(69) 



The driving force of the system seems to be equa- 
tion (68). Upon canceling common factors, this equation 
reads, from (53), 



-r(2 + a A + a gh - a/2) T{s/2 - q gh ) 
r(l + a gh ) T(s - a g h - a A - 1) 

= r(2 + a A + «q - a/2) Tjs/2 - q ) 
r(l + ao) r(s - a - o>a - 1) 



(76) 



We eliminate a A using a A = — 2a g \ L + (s — 4)/2 and 
obtain 



-r(-Q gh ) r(s/2 - a gh ) 

r(l + a gh ) r(( S + 2)/2 + « gh ) 

r(a - 2a g h) r(s/2 - a ) 



(77) 



T(l + a ) r((s + 2)/2 + 2a gh - a ) ' 



Note that 



-r(-agh) 



r(i - a gh ) 
a gh 



(78) 



9 



so this factor is positive for < a g h < 1. 
Now let s be fixed in the interval 

1< s < 3. (79) 

and let ao be fixed in the interval 

< a < s/2 < 3/2. (80) 

Then for a g h in the interval 

< a gh < a /2 < 3/4 (81) 

both sides of (77) are positive and finite. Moreover, by 
(78), when a g h approaches its lower limit, namely 0, the 
LHS of (77) approaches +oo while the RHS is finite, and 
when a g h approaches its upper limit, namely ao/2, the 
RHS approaches +oo while the LHS remains finite. Since 
there are no further poles or other discontinuities present 
in the stated interval, for every ao £ (0, s/2) there exists 
(at least) one solution a g h- 

This tells us that to solve (77) numerically we should 
take «o as the independent variable and we are assured 
that there exists a solution for 

<%h = a s h(ao), (82) 

both variables being in the stated intervals. 



which, with e = s/2 — ao, gives 

a g h(ao) = s/2 — a a « s/2 — e. (88) 

We have now determined that a g h(ao) vanishes at ao = 
0,s/2, and we have determined its slope at these two 
points. We approximate a g h(ao) in its interval by an 
interpolation. We write 

a gh = a (s/2 - a )/(a ), (89) 

where /(ao) is a positive function whose values at ao = 
and ao = s/2 are determined by the slope of a g h(ao) at 
these two points which are given respectively by 1 /3 and 
-1. 

A linear interpolation for /(ao) yields an approximate 
solution to (77), 

a gh (a ) «a (| - a ) + > ( 90 ) 

and comparison with the numerical results in subsec- 
tion VIII E shows that this is a reasonable approximation 
(with less than 10% deviation) even for intermediate val- 
ues of ao- 

E. Numerical Results 



D. Analytic Statements 

We may in fact solve (77) analytically for ao close to 
its end-points, ao = e and ao = s/2 — e, where e is 
small. Suppose first that ao = e. Then the inequality 
< a g h < e/2 implies that a g h is also small. In this 
limit (77) approaches 



-r(-a gh ) r(s/2) _ r(e - 2a gh ) T(s/2) 



r(i)r(( s + 2)/2) r(i)r(( s + 2)/2) 



(83) 



With — r(— a g h) « l/a g h, and T(e— 2a gh ) « l/(e— 2a gh ), 
we may equate the singular pole terms 



1 



1 



(84) 



a g h e - 2a g h 

and with e = ao, this has the solution 

a g h(a ) = a /3 a « 0. (85) 

Now suppose that ao = s/2 — e. Then since the RHS 
of (77) blows up like 1/e, a g h must be small, so (77) 
approaches 



-r(-a gh ) r( s /2) _ r( s /2) r( e ) 



r(i)r(( s + 2)/2) r(i + s /2)r(i) 

We again equate singular pole terms 

1 1 

a g h e 



(86) 



(87) 



A full analytic solution of (68) and (75) has not been 
obtained so far. To find a solution at least numerically, 
one can follow one of two strategies: 

• One equation (preferably (68)) can be solved for 
one variable (yielding a^ (ao)) and this function 
can be substituted into both F\ from (71) and F2 
from (74), yielding two functions of one variable 
-Fi(ao, a g h(ao)), i = 1,2. Intersection points of 
these functions are solutions of the system of equa- 
tions. 

• Equations (68) and (75) each implictly define a 

(possibly multi-valued) function a^{ao), i = 1,2. 
One can obtain each of these functions separately 
and find the solutions of the system as intersection 
points in a two-dimensional plot. 

The second strategy is more cumbersome, yet it gives 
a better understanding of the conditioning of the system 
and the relationship between solutions for various values 
of s. While a solution for (68) is easy to find, solving (75) 
requires more effort, since it defines a multivalued func- 
tion with several potentially relevant branches. 

We have employed the f indroot routine of Mathemat- 
ica 5.2 and 7.0.0 with a wide variety of initial guesses in 
order to find all branches. (In the color version differ- 
ent starting points can be recognized for having different 
shades of blue and green.) The corresponding plots are 
given in Figures 2 to 4. For comparison also the Fi vs. 
F2 plot for these three cases is given in Figure 5. 
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«0 


a gh 


QA 


a*/ 


s = 3 


1.07945 


0.240044 


-0.980087 


-0.599366 


s = 2 


0.351045 


0.105460 


-1.21092 


-0.140125 


s = 1 


0.256229 


0.068302 


-1.63660 


-0.119625 



TABLE I: Infrared exponents, obtained from the smallest so- 
lution of equations (68) and (75) for different values of s. Note 
that from (66) we have a m = a g h and a v = av . 





Q (1) 



(2) 


"o 


a (1) 


s = 1.06 






0.512478 




s = 1.05 


0.131834 


0.202188 


0.510570 


0.0408666 


s = 1.04 


0.093892 


0.224380 


0.508586 


0.0298297 


s = 1.03 


0.065958 


0.236660 


0.506537 


0.0212898 


s = 1.02 


0.041940 


0.245122 


0.504421 


0.0137068 


s = 1.01 


0.020197 


0.251402 


0.502242 


0.0066706 



TABLE II: Infrared exponents, obtained from the smallest 
of equations (68) and (75) for different values of s close to 
s = 1. The exponent a g h is only given for the solution close 
to Qo = 0. (Note that relation (85) is fulfilled quite well for 
this solution.) 



The two-dimensional plot is particularly interesting in 
the case s = 2 (Figure 3), since it shows that the sys- 
tem of equations is relatively ill-conditioned and the in- 
tersection point (i.e. the solution) would be sensitive 
even to small perturbations. (Such perturbations are of 
course absent in the present truncation, but could be in- 
troduced, for example, by vertex dressing.) 

For s = 3 and s = 2, the physical values are expected 
to be given by the smallest solution for ao, as summarized 
in Table I. 

For s = 1 the situation is somehow different since 
strictly speaking our equations are not well-defined (be- 
cause the transverse projectors vanish for s = 1). So 
instead of plainly taking the value obtained for s = 1 
we instead study the solutions for s = 1 + e with some 
small, but positive number e. For s = 1.06 there is only 
one solution at ao = 0.512478. 

When we lower s to s « 1.0541071910, in addition 
to the previous solution (which has now moved to ao ~ 
0.51136) another solution arises at ao « 0.17026. This 
solution splits into two independent solutions for smaller 
values of s; the smaller solution approaches zero for s — > 
1 + , as indicated by the values given in Table II. 

This is also illustrated in Figure 6 which, together 
with Figure 5.(c) gives a good impression of what is go- 
ing on: For s < 1.0541071910 there is a region around 
a Q k 0.17026 where F 2 > F x . For s = 1 + e with 
< e < 0.0541071910 we have F 2 (0) = and Fi(0) > 0, 
so there has to be an intersection point. For s = 1, 
however, we find F 2 (0) > F(0). The family of func- 
tions {F2} s= i+e. £ >o seems to converge pointwise, but not 
uniformly towards 1 s=i when e approaches zero from 
above. 



IX. RELATIONS AMONG THE 
^-COEFFICIENTS 

Having determined the 6 infrared critical exponents, 
the a's, we turn to the 6-coefficients. There are 6 DS 
equations and 6 coefficients 6a, b n i, 6 g h, bo, b m , b v . How- 
ever, as we have seen, only 4 of the equations for the 6's 
are independent, so there are 4 relations satisfied by the 
6's. These are, from (45), 



from (44) 



eq. (71) 



and eq. (69), 



&A&gh 



-1 



Is(V,S)i a A, «gh) ' 



1 



l V(V,S) 



(a a, a ) 



bo, 



F 2 (ao); 



(91) 



(92) 



(93) 



bl h = b b lp + b 2 m . (94) 
The last two equations determine the triple products 

1 



b A b Q b v 



F 2 {a)I V ( VtS - ) (a Al a<o) 



6 a tiL = 



-1 



F 2 (a)I V ( VtS )(a A i ao) 4(v,S)(«a, a g h) 



(95) 



(96) 



X. RANGE OF VALIDITY 

So far we have only taken into account the zeroth Mat- 
subara frequency, thus working effectively in the infinite- 
temperature limit. However, since infrared properties of 
the theory are governed by the zero frequency contribu- 
tion also at finite temperatures (for all T^0, even those 
in the confined phase) the results obtained so far may 
more general than initially stated. 

All Matsubara frequencies cu n with n > 1 effectively 
introduce mass terms in the non-instantaneous propaga- 
tors, so all such contributions decouple from the deep 
infrared where critical exponents are valid. Also in loop 
terms massive contributions show up only pairwise, so 
taking into account only the th Matsubara frequency 
yields a closed system. 

The only instance where zero and nonzero Matsubara 
frequencies could directly be intertwined is imposing the 
horizon condition - an issue which should be the subject 
of closer investigation. 
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FIG. 2: Solution of equations (68) and (75) for s = 3: We plot Q g h(o!o) from (68) and from (75). The relevant solution of the 
first first equation is represented by a solid black line. The second equation gives rise to a multivalued function: The graphs 
are composed of single dots; different shades of blue and green (in the color version) indicate different initial guesses. Note that 
single points which do not belong to any branch of the function typically correspond to values in which the f indroot routine 
got stuck. 

0f ell 

0.6r | 



0.4 



0.2 




FIG. 3: Solution of equations (68) and (75) for s — 2, otherwise as in Figure 2. 



XI. DISCUSSION AND SUMMARY 

A. Discussion of the Infrared Exponents 

We now turn to the discussion of the results given in 
table I. The horizon condition tells us that the ghost 
propagator D g h ~ 2 b fj a is enhanced in the infrared, 
so we are interested in solutions which fulfill 

ctgh > . 

As one sees from the black curves in Figs. 2, 3 and 4, this 
is true for our solution in the whole interval < a a < | 
which is the interval of physical interest. 

The most interesting quantity in Coulomb gauge Yang- 
Mills theory is presumably the color-Coulomb potential, 
given in momentum space by 

bnm 2a ° 

^Wfc)~7^Tw ( 97 ) 

It is linearly rising in position space for oio\ s = (s — l)/2. 

For s = 3, this gives ao|3 = 1, and the result we ob- 
tain lies slightly above this value, at = 1.07945. This 



corresponds to a slightly more than linearly rising po- 
tential 6 . The transverse gluon propagator -Daa(^) van- 
ishes at k = if a a < — 1. From Table I we have 
oa|s=3 = —0.980087, which corresponds to L>aa(&) 
that is weakly divergent at k = 0. This value is close 
to a a = — 1; which corresponds to a finite value for 
Daa(^ = 0), and a small change caused by an improved 
truncation could also change this to a a < — lj corre- 
sponding to a gluon propagator -Daa(^) that vanishes at 
k = 0. 

For s = 2, the value from Table I, = 0.351045, lies 
below the linearly rising case, = (s — 1)/2 = 0.5. How- 
ever as can be seen from Figure 3, the system is badly 



Note that the color-Coulomb potential is not a gauge-invariant 
quantity, so the arguments which forbid more-than linearly rising 
potentials in relativistic field theory do not directly apply; also 
the asymptotic inequality between Wilson and color-Coulomb 
potential in [19] would be satisfied. As discussed in section V, 
the mathematical aspects of such potentials are under control. 
Still it would be very surprising if in reality -Da j4 were more 
than linearly rising. 
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FIG. 4: Solution of equations (68) and (75) for s = 1, otherwise as in Figure 2. 



conditioned so that, if properly dressed vertices or terms 
neglected in our truncation even slightly modify the sys- 
tem of equations, then any value ao £ (0, 0.55) [which 
determines a value of a g h G (0, 0.146)] could qualify as a 
possible solution. This includes the linearly rising case. 
From Table I we find a a = -1.21092 < -1, which im- 
plies that Daa(^') vanishes at k = 0, but the uncertainty 
of ao, due to ill-conditioning, extends to all critical ex- 
ponents, including a a- 

For s = 1, the system is exactly solvable analytically. 
Since our equations do not strictly apply at s = 1 (see 
discussion in Sec. VIII E) we have examined numerically 
the limit s —> 1 + (sec Tabic II and Figure 6). The results 
obtained this way are consistent with eq. (85) which is 
valid for small oq. For s — > 1 + this solution converges 
towards 

a o = a gh = 0. (98) 

These values agree with the analytic solution presented 
in Appendix B. By (66) this would correspond to a^ = 
— |, but transverse propagators are not well-defined for 
s = 1, so this infrared critical exponent is undefined for 
the theory at s = 1. 

Because of the symmetry (47) the coefficient b is un- 
determined in the infrared asymptotic limit and must 
be fixed by subdominant terms that we have neglected. 
This is unfortunate because according to (38) and (39), 
for ao = 1, the quantity bo(g 2 T) 2a ° represents the color- 
Coulomb string tension and one would have liked to make 
a comparison with lattice determinations of this quan- 
tity [20-24]. 



B. Conclusion 

An ideal theory of confinement would explain in simple 
terms why there is a linearly rising Wilson potential and, 
in the Coulomb-gauge scenario, why there is a linearly 
rising color-Coulomb potential, with a Coulomb string 
tension a cou i (T) that increases with T even in the decon- 
fincd phase [20]. This we have not achieved; our solution 



of the Coulomb-gauge DSE at high temperature is ob- 
tained from equations involving r-functions that must 
be solved numerically. Nevertheless it remains true that 
the numerical value obtained for s = 3 space dimensions 
is numerically close to a linearly rising potential, with 

V(r) ~ r 2+2a -s = r l. 15890 

In s = 2 spatial dimensions the agreement is not as 
good but this may be due to the ill-conditioning of the 
equations. Finally, as s approaches 1 spatial dimension, 
the smallest solution approaches the exact analytic re- 
sult. 
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APPENDIX A: EVALUATION OF POWER-LAW 
INTEGRALS 

Determination of the infrared critical exponents re- 
quires evaluation of the integrals (48) to (51) which yields 
the results stated in (53) to (58). In order to illustrate 
the calculational scheme, we explicitly discuss the evalu- 
ation the integral (48). [The evaluation of such integrals 
is also discussed for example in the appendix of [25] ; here 
we give some more details on intermediate steps.] 

First we check issues of convergence: We write 

I s(v . s) (a A , <x v ) = h (A1) 

with 

j = f d s p k 2 p 2 -(p-k) 2 

1 ' J (271-)* (p2)2+a A [(jfc_ p )2]l+a„ " <■ > 
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(a) 



s = 1.05 




(b) 



1.01 



o.s 
0.6 
0.4 
0.2 / 



(b) 



s = 2 




(c) 



WFi) 
1.0 



1.002 



(d) 



s = 1.0004 



1.0 

o.s 






^1 


0.6 

0.4 








0.2 






1 


0.0 0.1 0.2 


0.3 0.4 


0. 



(c) 



S = 1 




FIG. 5: Solution of equations (68) and (75) for (a) s — 3, (b) 
s — 2 and (c) s = 1: We display i 7 ! (solid) and F2 (dashed) 
for as functions of ao with a g h(ao) from (68). We have em- 
ployed the mapping $(2;) = — arctan ^ to compactify the 
range (—00, 00) to (—1, 1). Solutions are determined by in- 
tersection points. Note that horizontal lines correspond to 
odd-order poles, so there is no solution for s = 3 at ao w 1.31 
and no solution for s = 2 at ao ~ 0.605. 



FIG. 6: Solution of equations (68) and (75) for (a) s = 1.05, 
(b) s = 1.01, (c) s = 1.002 and (d) s = 1.0004, otherwise 
as in Figure 5. One clearly notes that (as also suggested by 
Table II) there is a solution which approaches ao = a g h = 
for s — > 1 + but vanishes for s = 1. 



In this integral we obtain the power p s+2 from the numer- 
ator. In the infrared (p — > 0) the term (k~p) 2 is finite for 
k ^ 0, so we can neglect it for questions of convergence. 7 
Thus the denomiator contributes the power p 2 ( 2 + a ^) and 
a necessary condition for infrared convergence is 



s + 2 > 4 + 2a A 



(A3) 



(A2) may also have a non-integrable singularity at p = k. 
By the substitution p — > p + k the integral takes the form 

d s p k 2 p 2 -{p-kf 

(2tt) s [(p + jt) 2 ] 2 +«A (p2)i+a„ • \ > 

The denomiator now contributes the power p 2 ( 1+a v) an d 
we sec that we also have to demand 



s + 2 > 2 + 2a 



s 

<2 



(A5) 



7 Note that for k = the numerator vanishes, so this case is not 
problematic; for all other cases the argument holds. 
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in order to have infrared convergence. In the ultraviolet, 
k is negligible compared to p and the denomiator in (A2) 
contributes the power p 2 ( 2 + Q A+2(i+a ¥ ,)) _ -j nave ultra- 
violet convergence we have to demand 

g 

s + 2 < 6 + lax + 2a v , i.e. oa + Q v > - - 1. (A6) 

These conditions are compatible; by combining them one 
can also deduce the conditions a v > and a a > — 1. 

Now we proceed with the evaluation of the integrals. 
In (A2) we represent the powers of propagators by inte- 
grals over auxiliary variables. This employs the integral 
representation of the gamma function, 



k x J k x 
o 



> (A7) 



(fc > 0) which allows us to write 



1 



,1+CCA 



(p 2 )2+ a A J q da r(2 + a A ) 

1 f°° b a * 

db - 



-ap 



[{k-pYY+^ J T(l + a v ) 
With these identities 1\ takes the form 



(A8) 

b (k-p) 2 ( A9 ) 



It, = I da — / db „ 

/„ r(2 + a A )7o T(l + a v ) 



with 
h 



(2tt) s 



[fc V - (p • *) 



2 1 e"* 



/ 2 (A10) 



(AH) 



$ = ap 2 + 6(fc - pf = {a + b)p 2 + bk 2 - 2b(k ■ p) 

= (* + b ) (p-^ fc ) 2 ^ fc2 + &fc2 

= :( a + % 2 + -^-fc 2 . (A12) 

The substitution p — > g = p — A: yields 



#? (fcV-(<z-fc) 2 )c- (a+f 



(27T) 

-£^fc 2 ,,2 /" ds 9 



e"^ K fc' 



(2tt)< 



(l-(q-k) 2 )q 2 C^ a+b ^ 2 



One can rewrite an integral of the form 

J = k 2 f^-(l-(q-kf)f(<l 2 ) 



as 



(2ny 



J-fc.fc. I — 



(A13) 



(A14) 



-qiQj)f{q 2 



d s q 
(2n) s 



(Sij - /(g 2 ) = C ^ 



(A16) 



Taking the trace of this equation (6 a = s, q\q\ = 1) and 
a small rearrangement of factors yields 



C = 



s - 1 f d s q 2 
(2^ /(<?) 



So we have 



s - 1 



(A17) 



(A18) 



with 



h= l^o 2 e-^ 2 



(2tt) s 



1 



s+l -(a+b)q 2 



The substitution x = (a + b)q 2 yields 
1 



(A19) 



h = 



(2ny 
1 



dn s 

dtt s 



°° < x \ (s+1)/2 e"*dz 



K a + bJ 2x^2(0 + 6)V2 

1 



x 



(2tt) s 7 s 2{a + b) 1 + s / 2 J 

1 2 7r s / 2 1 fs 

rr - + 1 



(./a+D-ig-x^ 



(27T) 5 r(|) 2(a + 6) 1 + s / 2 V2 
1 1 f r (§) 



(4tt) s / 2 (a + fc)^/ 2 r(f) 
s 



2(47r)s/ 2 (a + fo)^/ 2 
Consequently we obtain (with slight rearrangements) 

(s-l)fc 2 e"^ fc2 



(A20) 



2(4tt) s /2 (a + 6)s/2+i 

To proceed with the evaluation of 

(s - 1) k 2 f°° , a 1+QA 
= — ; ; — t^- / da — ; 

2(4^/ 2 J Q r(2 + a. 



(A21) 



X / db 



a A 



e "+ b 



T{l + a v ) {a + byi 



2+1 



(A22) 



we introduce another auxiliary variable r by integrating 
over <5(a + 6 — r), 



(s - 1) fc 2 f 00 , f°° , a 1+Q 

■M = — ; ; — t^- / dr / da — ; 

2(4tt)-/ 2 Jo J q T(2 + c 



A 

a a) 



/o r(l + a^J (a + 6)-v 2 +i 



(A23) 



(A15) The substitution a = ra, b — t(3 yields 



Since J is a scalar which is quadratic in k, one has to i l 
find 



OO />OC 



(s-l)fc 2 1 

2(4^) s / 2 r(2 + a A )r(l + a v )7 " J 



(a + /J) 5 /^ 1 



(A24) 
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with 



/>OC 

Jo 



(A25) 



(Note that S(ra + t/3 - r) = i£(a + /3 - 1).) We now 
substitute t = a/3 k 2 r, 



(af3k 2 ) a * +a * +2 ~^ 7 
r(a A + a v + 2- f) 



and obtain 



( J^(aA+a p +2-f)-l e -t 

(A26) 



^s(v,s)(«A,a v ) 



iV(s-l) T(a A + + 2 - f ) 



2(4^/ 2 r(2 + a A )r(l + a v ) 



(A27) 



with 



CJC />oo 



(A28) 

Because of the delta functional, one only has contribu- 
tions from a + (3 < 1, accordingly one can restrict the 
integrations in (A28) to the interval [0, 1]. Performing 
one integration explicitly yields 



B V2 ~~ 2 ~ 1 ~ aA 
r(f -a v )TU -l-a A ) 



(A29) 



T(s — 1 — a v — Q!a) 
Plugging this result into (A27) yields 

I si v,s)^A,a gh ) = r(2 + aA) 

T(f -q y )r(f -l-a A ) 

r(i + a g h) r(s - 1 - a g h - «a) ' 

(A30) 

which is the result stated in (53). The other power-law 
integrals are evaluated in a similar way. Note that the 
integrals Iv(-,-) contain a factor (s — 1) _1 which stems 
from taking the trace of the projector on the lhs of the 
corresponding DS equation, 

trP«(p) = tr -^=5 u -^=s-l, (A31) 



and is included for convenience in the definition of the 
integral. 



APPENDIX B: ANALYTIC RESULT FOR s = 1 

We start from the second-order action (4) in s = 
space dimensions, 



1 



1 



5i = - J d Xl (jODiA,) 2 + d^cj . (Bl) 

where D\ (A) = d\ + gAi x , and the Coulomb gauge con- 
dition reads d±Ai = 0. Thus A\ is independent of X\. 
The propagators of A^ and the ghost are given by 



(Ao(xi)Ao(yi)) 
(c(xi)c(yi)} 



T 
T 



(B2) 



where the average is with respect to A\. Since A\ is 
independent of x\, we may diagonalize these operators 
by Fourier transform, so 



D Ao A {k\) = 
D gh (h) = 



(fci +i.gA 1 x) 2 
T 

kx(k! +igA 1 x] 



(B3) 



We now integrate A\ over the Gribov region, which is 
the region where the eigenvalues of the Faddeev-Popov 
operator k\(k\ +igAix) are non- negative. For this pur- 
pose, we quantize in a periodic box of length L, so 
k\ = 2Tm/L, where n is an integer. We first consider 
the gauge group to be the SU(2) group. One easily finds 
that the eigenvalues of igAi x are given by and ±|<?^4i | . 
The case k\ = is trivial. The Gribov horizon is de- 
termined by the first non-trivial zero eigenvalue. This 
occurs for n = ±1, at \gA\ = 2ir/L. Thus A\ is in- 
tegrated over the sphere \A\\ < 2-njgL. We now take 
the infinite- volume limit L — > 00, while keeping a typical 
momentum k± finite. In this case we have A\ — > 0, and 
the Aq and ghost propagators approach their free val- 
ues, Da Ao (fo) = £> g h(fci) = T/k\. The case of SU(iV) 
is the same. Accordingly the infrared critical exponents 
are given by «o = Qgh = 0- The restriction to the Gribov 
region (which in this case coincides with the Fundamen- 
tal Modular Region, where gauge-fixing is unique) was 
essential in deriving this result. 
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